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ABSTRACT 


For  large  values  of  the  pressure-gradient  parameter  0  (i.e., 
highly  accelerated  flows),  the  order  of  the  similarity  equations  gov¬ 
erning  the  laminar  boundary  layer  is  reduced  by  two.  Coles  has  shown 
that  this  is  a  singular  perturbation  problem  for  which  an  inner  solu¬ 
tion  and  an  outer  solution  must  be  obtained.  In  this  paper,  inner  and 
outer  solutions  are  presented  which  include  the  effects  of  a  power- law 
temperature- viscosity  relation,  nonunit  Prandtl  number,  leading-edge 
sweep,  and  surface  mass  transfer.  An  extremely  accurate  method  is 
described  for  estimating  skin  friction  and  heat  transfer  for  all  posi¬ 
tive  values  of  0. 


Any  views  expressed  in  this  paper  a*e  those  of 
should  not  be  interpreted  as  reflecting  the  views  of 
tion  or  the  official  opinion  or  policy  of  any  of  its 
private  research  sponsors.  Papers  are  reproduced  by 
tion  as  a  courtesy  to  members  of  its  staff. 
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SYMBOLS 


a  =  constant  of  the  inner  solutions,  o  cos2A/t 

2  3 

c  =  fluid  specific  heat  (taken  to  be  constant  across  the 
boundary  layer) 

f  =  variable  related  to  streamwise  velocity;  f  '  =  u/u 

e 

fw  =  surface  mass  transfer  parameter,  -  [(pv  )/(p  li  u  )  1  /TZ 

w  w  w  ey  J  ^ 

g  =  transverse  velocity  function,  w/w 
H  =  total  enthalpy,  h  +  u  /2 
h  *■  static  enthalpy,  c^T 

*1  (1  '  *  I1<2)]  *  (1  -  tw)Ii(3)  +  I  (2) 

00 

12  =  /f/<L  -  f/)dTl 

0 

00 

13  -  Jf'(l  -  f/2)dTi 

0 

oo 

^(1)  -  /(l  -  g2)dTj 

0 

1,(2)  -  |(1  .  f '2)dT| 

OO 

I]_(3)  =  /(I  -  8)dT) 

0 


k  =  thermal  conductivity 
M  ••  Mach  number 
Pt  °  Prandtl  number 
P  ■  pressure 
R  ■  gas  constant 

s0  "  transformed  velocity  function  defined  by  Eq.  (24) 
T  =  fluid  temperature 
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free-strear  stagnation  temperature,  H£/Cp 

transformed  similarity  variable  defined  by  Eq.  (24) 

sweep  parameter,  [l  +  { ( Y  -  1)/2}M^  cos2  A]/[l  +  {(V  -  1)/2}M(jd1 

wall  temperature  ratio,  Hw/H£  or  T^/Tq 

free-stream  velocity 

flow  velocity  in  the  x-direction 

flow  velocity  in  the  y-direction 

flow  velocity  in  the  z-direction 

constant  defined  by  Eq.  (15)  ’ 

coordinate  in  the  direction  of  flow 

coordinate  normal  to  the  surface 

coordinate  transverse  to  the  flow  direction 

pressure-gradient  parameter  defined  by  the  similarity  re¬ 
lation 


u 

e 


r  bi1/2 

L<Te/T0)5  J  ’  i'e-  0 


u 

e 


specific  heat  ratio 

similarity  variable  defined  by  Eq.  (3) 

dimensionless  total  enthalpy  (H  -  Hw)/(H£  -  H^) 

angle  of  sweep  between  leading  edge  and  the  normal  to 
the  free-stream  direction 

densi ty-viscosity  ratio  (pn/p^^) 

viscosity 

transformed  x-coordinate  defined  by  Eq.  (4) 


* 


l 


fluid  density 
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2 

a  ■  hypersonic  parameter  (U  / 2H  ) 

®  6 

2  9 

ct.  =  modified  hypersonic  parameter  (U  / 2H  )  •  (u  /u  ) 
i  ®  e  e  ® 

2  2  2  2 

ct9  -  modified  hypersonic  parameter  (U  /2H  )[(u  /u  )  cos  A  +  sin  A] 

4  ®  0  0  CD 

u>  *  exponent  in  the  temperature-viscosity  law  ^  ~ 

Subscripts 

(  )  *  value  at  the  edge  of  the  boundary  layer 

(  )  -  value  at  the  wall 

w 

(  )„,  ■  value  of  the  function  in  the  free  stream 
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I.  INTRODUCTION 

The  similarity  equations  governing  the  laminar  boundary  layer  are 
well  known.  For  a  thermally  and  calorically  perfect  gas  of  homogeneous 
composition,  there  are  three: 

St Teamwise  Momentum  Equation 


<Xf")  '  +  ff* 


p{f  ,2  *•  r  L(1  •  V6  •  (1  -  t8)g2  +  cw] )  <l> 


Transverse  Momentum  Equation 


(Ag ;) '  +  fg  '  -  0 


Energy  Equation 


e  y  +  te  * 


{(Htt  (g  -  0 


cos  A  +  gg  '  sin 


In  these  equations,  the  quantities  f',  0,  and  g  are  the  normalized  stream- 
wise  velocity,  total  enthalpy,  and  transverse  velocity  respectively,  and 
primes  denote  differentiation  with  respect  to  the  transformed  variable 
T|  normal  to  the  surface; 


T1 


u  y 

ml"’ 


where 


(3) 


/  ‘W‘e<lx 


0 


(A) 
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Other  quantities  are  defined  in  rhe  list  of  Symbols. 

These  equations  are  spplieable  only  when  the  boundary  layer  is  self, 
similar ,  i.e..  when  derivatives  of  f,  0,  a„d  s  with  respeet  to  S  are  aero. 
Often,  however,  it  is  possible  to  obtain  reasonably  accurate  values  for 
boundary. layer  properties  (including  the  shin  friction,  best- transfer 
rate,  and  displacement  thichness,  by  assuming  that  these  properties  are 
Slven  by  the  similar  solution  corresponding  to  the  local  values  of  0, 
tW’  ts>  etc- 

In  using  the  concept  of  local  similarity  in  highly  accelerated  flows, 

useful  t,  have  solutions  of  the  laminar-boundary- layer  equations 

for  large  values  of  0.  Thc  mathematical  difficulties 

an ilcul ties  encountered  in  the 

limit  3  »  i  may  be  inuatrated  .  a  . 

y  ng  the  streamwise  momentum  equa¬ 

tion  in  the  form 


5  W*V  +  if]  -  {f'2  -  i  [<i  .  tu)e  .  (I  .  t<)g2  +  g) 


o  (5) 


In  the  limit  3 

.  tq.  (5)  reduces  from  third  to  first  order: 

0~  £  ^  10  *  C“)9  •  (1  -  's^2  +  fwlj1^  (6) 

The  transverse  momentum  equation  and  the  energy  equation  remain  „f 

ond  order.  As  originally  pointed  out  by  Coles<2>  and  as  later  elaborated 

by  Beckwith  and  Cohen  thli  in,Hs  a 

n,  this  leads  to  a  singular  perturbation  problem 

in  which  the  thickness  of  the  velocity  layer  is  of  order  0‘ 1/2  „lth 
-Pectoral  enthalpy  a„d  transverse  velocity  layers  of  order  unity.* 

Discussion  of  this  problem  also  appears  in  an  abbreviated  version 
in  Lagerstrom's  article.'4'  version 
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The  appropriate  method  of  .elution  1.  si,ll.r  to  that  employed  in  de¬ 
riving  a  uniformly  valid  approximation  to  the  Navier-Stokes  equation, 
m  the  limit  of  large  Reynolds  number  (see  Kaplun  and  Ugerstrom, (5) 
Lagerstrom  and  Colc,<6>  and  the  reeent  book  by  Van  Dyke(7)). 


The  Outer  Limit  Equations 

Since  the  mathematical  Justification  of  this  singular  perturbs 


tlon  solution  (more  popularly  called  an  Inner-  and  outer-expansion  pro¬ 
cedure)  has  been  discussed  In  detail  by  Coles, <2>  our  purposes  will  be 
served  by  a  cursory  development  of  the  governing  equations.  The  pre¬ 
sent  analysis  extend,  the  work  of  geckwlth  and  Cohcn(3)  to  include  a 
Power-law  temperature-viscosity  relation,  a  constant  but  nonunit  Prandtl 


number,  and  surface-mass  transfer. 

Assume  that  appropriate  outer  representation,  of  the  dependent 
variables  f,  0,  g  are  Qf  the  following  fonng; 


f  “  fo  fi  +  ••••! 

6  "  eo  +7b  ei  + 


8  "  8o+/p-gi  + 


(7) 
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Substituting  these  representations  into  the  momentum  and  energy  equa- 
tions  and  dropping  all  terms  of  order  8' 1/2  and  smaller,*  the  following 
"outer"  equations  are  obtained: 

£o'  ■  if  {<*  -  ‘w>80  -  (1  -  +  C„}J1/2  (8) 

s 

(ir  eS)'  +  f  e'  +  ( — Ci  - 

VPr  0 J  0  0  1(1  -  t  )  V  Pry) 

w 

~2~  (JQ  COs2  A  +  ~2  Sln2  A]  }  “  0  (9) 

'  +  f080  *  0  (10) 

The  appropriate  boundary  conditions  are  found  ^a)  by  requiring  that 
the  outer  equations  satisfy  the  outer  exact  boundary  conditions,  and  (b) 
by  exact,  matching  of  the  inner  and  outer  representations.  The  results 
may  be  written 

V°>  “  V  60<°>  "  80<0)  =  °»  0o<®)  "  *  1  (H) 

Only  one  boundary  condition  on  fQ  may  be  satisfied  by  the  outer  equa¬ 
tions,  because  the  outer  limit  equation  for  fQ  reduces  to  first  order.  The 
exact  boundary  condition  f  -*  1  as  T)  -♦  oo  is  automatically  satisfied  by  Eq. 
(8).  The  no-slip  condition,  f '  -  0  as  T|  -•  0,  must  be  satisfied  by  an 

* 

The  outer  limit  equations  are  properly  obtained  by  applying  the 
limit  8  -*  00  to  the  full  equations  expressed  in  outer  variables,  with 
1  held  fixed.  This  gives  identical  results  to  those  cited  here. 


inner  solution  which  is  valid  in  a  region  of  extent  3" 1/2  with  respect 
to  the  scale  of  the  outer  solution. 

In  this  approximation,  the  density-viscosity  ratio  X  ■  (p^/p^pi  ) 
in  the  outer  layer  becomes 


1  -  CT, 


l)  -  Veo  - 


U)-l 


(1  - 


+  t  } 
wJ 


(12) 


where 


and  we  have  assumed  that  the  viscosity  is  equal  to  a  power  function  of 
the  temperature,  i.e.  ,  p  ~  I*.  M-I,N[(i.  „,)/t  because 

f0(°)  *  and  f0(°)  ls  nonzero  in  general.*  The  energy  and 

transverse  momentum  equations  are  coupled  through  the  implicit  appear¬ 
ance  of  both  eQ  and  gQ  in  Eqs .  (9)  and  (10).  If  Pr  -  1,  eQ  -  ^  and 
the  number  of  coupled  ordinary  differential  equations  is  reduced  from 
three  to  two.  Setting  both  the  viscosity- temperature  exponent  uu  and 
the  Prandtl  number  Pr  equal  to  unity  reproduces  the  equations  of  Beck¬ 
with  and  Cohen. ^  ^  If  -  0  (i.e.,  the  local  Mach  number  is  zero), 

ts  -  1  and  the  transverse  momentum  and  energy  equations  are  again  un¬ 
coupled. 


It  should  be  noted  that  the  inner  and  outer  expansion  procedure 
breaks  down  in  the  limit  ty  -•  0,  because  the  outer  solution  for  fn 
satisifies  the  exact  boundary  conditions  for  the  complete  equations 
and  the  inner  solution  for  f'  is  simply  zero. 
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One  interesting  case  was  pointed  out  by  Coles  and  we  extend  his 

result  to  generalized  compressible  flow.  For  t  »  t  ■  1.  A  -  0  so 

w  s 

that  a2  -  (uf/2He)(ue/uJ2  e  Oj,  X  -  ( 1  -  a^"1,  and  fQ'  -  1  for  all 
T\.  In  taking  this  limit  with  Pr  +  1,  the  product  of  (fQ/2)  #  and 
approaches  9^  so  that  the  energy  equation  becomes 

eS  +  (71  +  *w>Ml  -  -  Pr)]-1  (l  -  a/"®  9(5  -  o  (i3) 

with  the  boundary  conditions 

eQ(0)  -  0;  9q(oo)  -  1  (14) 

If  we  define  the  new  variable  x  and  the  constants  Xq  and  W  by 

V  m  —— 1  _  V 

x  rt  xo' 

*0  (15) 

(i  - 

w  "  U  -  aL(l  -  Pr)]Pr 

then  Eqs.  (13)  and  (14)  are  satisfied  by  the  solution 

9  -  Cerf(X//I)  +  erf  (x^/I)  ]/[l  +  erf  (x^/?)]  (16) 

For  Xq  <  0,  we  note  the  identity 

erf  (-x)  -  -erf  (x)  (17) 

The  corresponding  solution  for  gQ  is  found  from  Eq.  (16)  by  setting 
Pr  -  1. 
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Equations  (8)  to  (10)  with  the  boundary  conditions  of  Eq.  (11) 
represent  a  two-point  boundary. value  problem  for  the  complete  solution 
of  the  outer  limit  equations.  The  derivatives  6^(0)  and  g^O)  ,  along 
with  the  Integrals  I,.  1,(1),  1,(2),  and  1,(3)  evaluated  using  f  « 
and  gQ  In  place  of  f,  0,  and  g,  are  given  In  Table  1.  The  Integral  I 
Is  Identically  sero.  Beckwith  and  Cohen<3>  calculated  several  „f  thel 

quantities  for  the  special  case  of  Pr  .  <e  .  1,  f^  .  0,  and  nonunit  val- 

ues  of  t  . 
s 

The  Inner  Limit  Equations 

Inasmuch  as  the  order  of  the  energy  and  transverse  momentum  equa¬ 
tions  are  not  reduced  in  taking  the  limit  0  -  »,  the  outer  equations 

(Eqs.  (9)  to  (11)),  represent  complete  solutions  for  the  total  enthalpy 
and  transverse-velocity  profiles  for  large  0.  The  terms  fQ  and  (fj2) , 
which  appear  in  the  outer  equations,  differ  from  the  exact  solutions  f 
and  (f/2)  only  in  a  region  which  is  0'1/2  smaller  in  extent  than  the 

region  of  applicability  of  the  outer  equations.  Therefore,  the  outer 

equations  asymptotically  represent  the  complete  solutions  for  0  and  g 

as  0  -  o»,  and  ln  thl8  llmlt  the  lnner  8olutions  foj.  Q  flnd  g  aj_e  ldentl_ 
cally  zero. 

The  no-sllp  condition  f '(0)  -  0  1.  setlsflcd  by  the  Inner  llmlt 

equations  for  f;  the  Inner  equation,  for  0  and  g,  as  noted  previously, 

reduce  to  0  .  g  .  0.  To  examine  the  Inner  streasvl.e  nxmentum  equation, 

It  Is  necessary  to  Introduce  a  new  Independent  variable  if  and  an  Inner 
representation  for  f: 


tT  -  yrn;  f 


+  . . . 


(18) 
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so  Chat 


(19) 


After  rearrangement,  the  inner  equation  for  f '  becomes 


(TJ>2  +  r 

s 


(20) 


and  dropping  terms  of  order  1//F  and  smaller,  the  result  is 


where 


(*?0)  ' 


<*o>  +  r 

s 


X 


[i 


2 

cos  A 


t 

v 


(21) 


(22) 


The  boundary  conditions  on  fQ'  are  found  from  the  exact  boundary  condi¬ 
tions  at  the  wall  and  the  matching  condition*  that  the  inner  and  outer 
representations  of  f  agree  in  the  limit -  “  with  if  large  but  fixed. 
The  following  boundary  conditions  for  the  inner  equations  result: 

f0(°)  "  0;  f0<°>  "  °»  V”*  “  V  ty/tg  (23) 


* 

The  matching  condition  applied  here  is  elaborated  by  Van  Dyke-^ 
Inner  representation  of  (outer  representation)  -  outer  representation 
of  (inner  representation). 
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Si  nc  e  Eqs .  (21)  and  (22)  involve  only  £q  and  not  f^,  they  repre¬ 
sent  a  well-posed  second-order  two-point  boundary- value  problem.  A 
more  symmetrical  form  may  be  obtained  by  applying  the  transformation 

80(t)  "  (tw/t:s)1/4fl  (24) 

Then  the  problem  may  be  written 

(Xfl0)  7  ■  8o  *  1  ’  0  (25) 

x  "  ll  -  <T  (26) 


s0(0)  -  0;  sQ(“)  -  1 


(27) 


where  the  constant  a  is  given  by 

2 

a  -  (a.  cos  A)/t 

i  s 


(28) 


Values  of  sQ'(0)  satisfying  Eos.  (25)  to  (27)  are  given  in  Table  2 
for  several  values  of  u>  and  a.  By  using  a  simple  transformation  sug¬ 
gested  by  J.  Aroesty  the  two-point  boundary- value  problem  of  Eqs.  (25) 
to  (27)  is  reduced  to  the  numerically  simpler  problem  of  a  double  quad¬ 
rature  from  0  to  t.  Define  a  new  variable  a  by  the  relation 


80  "  d80/dt  =  * 


(29) 


so  that  the  transformation  (sQ,t)  -•  (s0><*)  is 


The  RAND  Corporation,  private  communication. 
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_d_ 

dt 


(30) 


and  hence 


“  <^>  ■  4 ♦ 1 


(31) 


A  subsequent  transformation  (s^.a)  -*  (s^.T)  is  defined  by 


T  ■  \ct 


(32) 


and  upon  multiplying  Eq.  (31)  by  X,  we  have 

-dT  +  \(s0)[l  -  SQ]ds0  -  0  (33) 

with  the  single  boundary  condition  T(Sq  =  1)  *  0.  Separating  variables 

and  integrating  Eq.  (33)  from  s^  *  1  to  Sq,  we  have 


(34) 


A  final  quadrature  of  Eq.  (34)  suffices  to  compute  t(s^)  ;  noting  that 
T  »  X(dsQ/dt)  and  sQ(t  ■  0)  ■  0,  we  obtain 

t(s0)  -  J*  [*(.*)/T(.*)]da*  (35) 

0 

The  quadratures  of  Eqs.  (34)  and  (35)  are  rapidly  accomplished 
using  standard  numerical  integration  formulas.  The  skin- friction 
derivative  Sq(0)  is  found  from  the  quadrature  of  Eq.  (34);  because 
X(0)  =  0  for  all  uu  and  a, 
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8q(0)  =  a(0)  =  t(s0  =  0)  (36) 

For  the  special  case  of  to  =  1,  Coles^  ^  pointed  out  that  an  analytic 
solution  of  Eq.  (25)  may  be  found  in  the  form 

sQ  =  1-3  sech  (t//7  +  tanh"  V2/3)  (37) 

Using  either  the  analytic  solution  for  to  ■  1  or  the  numerical  solutions 
for  to  4  1,  the  surface  skin  friction  derivation  f*  is  then  found  by  re- 
versing  the  previous  transformations: 


lim  * 

g-MO  w 


(38) 


In  general,  8^(0)  depends 
to  ■  1,  the  value  of  8^(0) 


on  the  three  parameters, 
is  4/3. 


’2’ 


to,  and  t  ;  for 
s 
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■II.  discussion 

It  is  very  difficult  to  obtain  exact  numerical  solutions  of  the 
laminar  boundary- layer  equations  for  values  of  0  greater  than  2.  The 
reason  is  simply  that  the  singular  behavior  of  f /(T|)  near  the  wall 
which  exists  for  P  -*  ®  becomes  dominant  even  for  moderate  values  of 
3.  Conversely,  the  results  obtained  for  P  -  ®  should  be  good  repre¬ 
sentations  of  the  behavior  of  the  boundary  layer  for  large  but  finite 
values  of  p.  One  of  the  important  results  that  we  wish  to  demonstrate 
is  that,  by  combining  the  exact  numerical  results  of  Ref.  1  for  P  55 
and  the  limiting  solutions  for  P  -  «  given  in  Tables  1  and  2,  it  is 
possible  to  estimate  accurately  the  skin-friction  and  heat-transfer 
derivatives  f",  gj,  and  for  all  positive  values  of  P. 

The  skin-friction  results  for  P  ■  are  displayed  in  Fig.  1.  The 
influences  of  the  two  parameters  uo  and  a  on  the  inner  solutions  for 
•qCO)  are  seen  to  be  relatively  small.  For  0.5  ^  u)  <:  1.0  and 
0  <  a  <  1.0,  the  value  of  8^(0)  may  be  found  by  interpolation  to 
better  than  0.25  percent. 

The  skin  friction  parameter  f^P  ^(t^/t^)  approaches  the 
limit  of  Sq(0)  as  P  -*  ®.  The  difficulties  that  were  observed  pre¬ 
viously^  in  calculating  exact  solutions  for  p  £  2  imply  that  this 
limit  is  approached  very  rapidly  with  increasing  p.  This  supposition 
is  borne  out  by  Figs.  2  to  A,  where  the  skin  friction  parameter  is 
shown  as  a  function  of  p"1/2.  The  limit  parameter  p_1/2  is  suggested 
by  whe  ordering  procedure  used  to  obtain  the  inner  and  outer  equations. 
Solid  lines  indicate  exact  numerical  solutions  from  Ref.  1  and  dashed 
lines  are  extrapolations. 


Figure  2  illustrates  the  approach  of  the  skin-friction  parameter 

to  its  asymptotic  limit  Sq(0)  for  different  wall  temperatures.  The 

limiting  value  is  approached  most  rapidly  for  high  wall  temperatures. 

This  result  is  to  be  expected  from  the  behavior  of  the  outer  equations. 

Large  values  of  tw  increase  the  magnitude  of  the  velocity  difference 

across  the  inner  layer,  whereas  for  ty  -  0,  the  distinction  between 

the  inner  and  outer  layers  breaks  down  and  no  proper  limit  is  obtained. 

It  has  also  been  found  empirically^  that  exact  solutions  are  more 

difficult  to  obtain  for  large  t  . 

w 

The  approach  of  the  skin-friction  parameter  to  its  limiting  value 
with  increasing  3  is  illustrated  in  Fig.  3  for  several  values  of  the 
sweep  parameter  tfi.  From  a  numerical  point  of  view,  the  accuracy  of 
the  present  extrapolation  procedure  increases  with  increasing  sweep. 

It  is  of  interest  to  examine  the  behavior  of  the  inner  and  outer 
equations  with  mass  injection  at  the  wall.  Applications  of  these  re¬ 
sults  include  mass  transfer  cooling  of  rocket  nozzles  (Back  and  Witte 

(8)v 

)  and  blunt  hypervelocity  vehicles.  The  outer  equations  determine 
the  heat  transfer  derivative  6j  and  they  contain  the  boundary  condi- 
tion  fQ(T|  -  0)  -  fw»  Thus,  surface  mass  transfer  (f  <  0)  acts  to  re¬ 
duce  ej  and  consequently  surface  heat  transfer  even  in  the  limit  3  - 
On  the  contrary,  the  skin  friction  derivative  f'  is  found  from  the  in¬ 
ner  solution  for  f'(0),  Eq.  (21),  and  the  solution  of  this  equation 
is  independent  of  the  value  of  f^  In  highly  accelerated  flows,  there¬ 
fore,  the  effects  of  blowing  on  skin  friction  become  negligible  in  the 
limit  3  -  ®  with  f^  fixed. 

Figure  4  illustrates  the  behavior  of  the  skin-friction  parameter 
as  a  function  of  the  injection  parameter  f  . 


The  limit  is  approached 


smoothly  with  decreasing  values  of  P 


for  all  f  ,  and  accurate  es- 
w 

timations  of  f  may  be  obtained  for  all  P  by  comparing  the  exact  solu¬ 
tions  of  Ref.  1  for  P  i  5  and  the  inner  limit  solutions  of  Fig.  1. 

The  wall  heat  transfer  is  related  by  the  modified  Stewartson  and 

Howarth-Dorodnitsyn  transformations  (Eqs.  (3)  and  (4))  to  0  '.  Inas- 

w 

much  as  0^  and  g^  represent  the  complete  solutions  for  0  and  g  as 
P  the  values  of  0^(0)  and  g^O)  represent  the  asymptotic  limits 

of  0j  and  as  P  —  <*>.  A  comparison  between  our  present  results  and 
Ref.  1  suggests  that  the  asymptotic  values  are  within  a  few  percent  of 
the  exact  answers  for  P  ^  7.  It  should  be  emphasized  that  0j  and  gj 
become  independent  of  P  as  P  -*  «,  demonstrating  that  the  heat  trans¬ 
fer  predicted  by  a  local  similarity  analysis  is  highly  accelerated 
flows  approaches  a  limiting  value. 

Figure  5  shows  the  typical  behavior  of  the  heat-transfer  deriva¬ 
tive  9j  with  increasing  values  of  p.  The  behavior  of  0j  with  decreas- 
-1/2 

ing  values  of  P  is  seen  to  be  smooth  and  (at  least  for  the  case  of 

(9) 

■  0)  monotonic.  In  an  earlier  paper,  we  demonstrated  that  the 
proper  parameter  to  use  in  comparing  different  heat-transfer  calcula¬ 
tions  is  0/[(l  -  t  )/(t  -  t  )1.  For  a,  ■  0  and  t  ■  1,  the  adiabatic 

wall  temperature  t  is  unity  for  all  Pr  so  that  the  heat  transfer  pa¬ 
rameter  reduces  to  9j.  Although  we  have  not  been  able  to  prove  it  an¬ 
alytically,  it  appears  that  t  ■  1.0  for  all  values  of  Pr,  t  ,  uj,  and 

flW  8 

a  in  the  limit  P  -  ®,  This  is  a  very  surprising  result  and  should  be 
examined  further. 

Hr 

The  boundary-layer  displacement  thickness  6  is  defined  by  the  re- 
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(39) 


where  the  integrals  I  and  I  are  given  in  the  list  of  Symbols.  For 
0  -  ®,  the  integral  1^  is  identically  zero  and  the  integral  I 


I 


2 


f'(l  -  f')  dT| 

0 


(40) 


is  recorded  in  Table  2.  In  the  absence  of  sweep,  the  velocity  profile 
is  monotonic  and  0  <  f  '  <  1,  so  that  the  integral  I  is  positive  and 
the  displacement  thickness  is  negative.  With  sweep  (tg  <  1),  there 
is  often  an  overshoot  in  the  velocity  profile  so  that  f '  >  1  f or  some 
range  of  T)  and  I  becomes  negative.  With  sweep,  therefore,  the  dis¬ 
placement  thickness  may  be  either  positive  or  negative,  depending  upon 
the  particular  parameters  being  considered.  Numerical  comparison 
between  the  calculated  results  of  Ref.  1  for  moderate  0  and  the  pre- 
sent  results  for  0  «  »  suggests  that  6  monotonically  decreases  with 

increasing  0.  In  the  special  case  when  0  «  »  and  t  ■  t  *1  f  '  ■  i 

w  s  ' 

*  -1/2 

and  6  is  of  the  order  of  0  times  the  scale  of  boundary  layer  dis¬ 
placement  thickness  for  (3  «  0. 


SOLUTION  OF  THE  INNER  LIMIT 
EQUATIONS  FOR  g  » 
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SOLUTION  OF  THE  OUTER  LIMIT  EQUATIONS 


ml  sfss  S253  SlSs  SSSS  sSSs 


H  H  H  N  H  H  H  N 


2  n  oo  N 

°  9  o  —i  r*i 

O'  ®  w  o  <r  cn 

i— i  oo  in  o  o\  vo  4  o 


3  £  g  _  <r  <r  vooo 

r?  oi  2  o^n  '+2^  «  ^  O 

°1  ^  o  O  2  O  hS°o 


^  O  O  O  "-I  O  o  o 


n  «  o>n 
>a  ®  o  o 
oo  O'  o  o 


^  N  n  <-l  r-l  i— ( 


3£S;j  2  £  £  g  2S^°°  o^rvoo 

g  g  g  3  oh22  33SS  22S£ 


CM  CM  CN  cn 


P"*  VO  so 

oo  in  cm 

in  oo  r^ 

<n  cn  o 
•  •  •  • 

o  o  o  o 


noS  £  2  o  £  )2  2  O'  <f  °o 

<f  oo  m  i^o-h  osZiS  £0-3  h  ^  h 

-3v3cmo  SSSo  3££o 

oooo  OOOO  OOOO  OOOO  OOOO 


<}  -JN 
n  o'  \o  ^ 
f'  -j  to  n 
m  nn  <f 


VO  H  r-l 

O  IA  r-l  r-l 

n  >J  o  ffl 
n  -J  in  m 


g'OOn 

00  CM  -3  v£> 

o  H  ID  fl 

n  <j  vf  m 


oooo  oooo 


o  to  oi  m 

io  h  a  oi 

m  h  sj  s 

CM  CA  (A  O 


i-1  -J  lO  i-l 
n  id  h  o 

m  a  n  o 
n  n  ^  in 


— i  m  oi 

•J  a  a  oo 
-1  -3  oi  a 
m  n  n  sj 


o  o  o  oooo 


CM  00  O'  r- 
0\  kD  <t  OO 

n  in  O'  >j 
<n  on  m  <t 


VO  O'  r-l  r-l 

'O  CM  o  in 
r-i  oo 

tn  a  a  a 


<n  tn  a  a 
•~<  M3  on  on 
on  -o  o  oo 
n  •}  in  m 


OOOO  OOOO 


on  vo  un  cm 
O  a*  in  o' 
iJnino 
<n  on  on  <t 


O  N  O'  O' 
on  on  oo  _i 
n  o  a  a 

CM  -3  -3  IA 


OOOO  OOOO 


n  in  o  o' 
a  rv  n  -I 
IA  vO  CM  O 
CM  tn  <f  IA 


oooo  oooo 


r-l  -3  VO  o  r-l>3vOO 


r-l  -3  a  o 


m 

H>JiOO 


m 

H  ^  \D  o 


in 

-«  ^  v*>  O 


>*l3e 


Table  2,  cont 


-23- 


m  O'  n  r-'  os  cn 

<r  o  m  <r  m  *j  m  co  <j-  ri 

oo  oo  m  o  n  o  oo  n  o  oo 

so  so  so  so  n»vom  so  so  m 

••••  •••  •  •  • 

4  i-h i  eg  O  O  O  OOO 


on  <r  o 

r-  ct>  o 
vD  vt  <NI 
vO  vD  vO 


00  00 
moo 
O'  m  \o 
00  o  oo 


vO  O'  N 
r-l  <SI  O 
^  ON  N 
00  N  N 

•  •  • 
OOO 


N  OO  N 
O'  in  <j 
m  ^ 
oo  oo  oo 


vo  o  m 
'O  o  in  <y 
O'  o  oo  m 

O'  o  O'  O' 
•  •  •  • 
O  ^  O  H 


CT> 

H  OOh 

m  o  vo 
<r  o  v.0  o 

•  •  •  • 
H  r-M  O  O 


m  o  m 

O  OOlfl 
in  n  w 
•  •  • 
eg  m  <r 

i  i  i 


h  o  m 
n  oo  no 
eg  oo 
•  •  • 
eg  m  <t 

i  i  i 


no  vo  in 
Mn  n 
m  n  4 
•  •  • 
eg  <j-  m 

i  i  i 


vD  ^ 

O'  <f  m 
m  m 

h  co  n 
•  •  • 
OO—I 
I  I  I 


so  e'¬ 
en  o  in 
O  cm  r- 
O  SO  o 


o  m 
cm  o 
O'  >j  n 
O  I-"  CM 

•  •  • 
OO-i 
I  I  I 


— I  O  — I 

n-  o  <r 

£  o  os 

°o  so  cn  o 

*>ii 

o  o  o  o 


O'  n  co  n-rnoo  as  oo  in 
cm  in  so  n-  cm  m  as  n  O'  so  m  cn  as 

o  m  <f  o  f'  cm  as  vi  cm  as  h  o  to  o 

CM  CM  CM  CM  OS  CO  N  CO  CO  fN  Os  O'  00  CM 


O'  H  O' 
O  — c  00 
CM  CM  —I 


00 

'0  4  0'  as  so  -j  o  m  oo  oo 

n  N  O'  CM  CM  O  00  00  sO  CM 

O  O  O'  —I  — i 


CMCMCMcn  OOO  n  n  n  (— s  r-1  f— 1  _ ,  ,  ,  *  1  *  *  •  •  •  • 

^  c->oa  ^  — 1  — '  O  — i  — i  — i  — i  — c  — c  cm 


— <  Os  so 

in  as  so 
<r  cn  o 
m  4  rs  o 
•  •  •  • 
o  o  o  o 


_  CM  O  SO  H  00  00 

^2  ^  ^  ^  ^  cm  O'  o  n  4 

rj  *2  ^  -2  ^  tMfMvo  oo  m  m  -t  m  o 

^  -I  n  n  4  cm  mfMso  o  sj  ts  0  cn  so 


•— <  eg  rn 

i  i  i 


•— i  eg  m 

i  i  i 


*-*  oo  m 

m  ^  co  h  n  n 

m  o  av  n  no  oo 

o  <r  vo  n  cj  h  o 


ooo  ooo 

. 


ooo 

I  I  I 


o  o  o  o 


egooeg  co  eg  eg  av  eg  eg  eg  oo  vo 

^vpmm  vo  ^  oo  av^co  in  av  oo 

or^^vo  <M^^h  oo  <f  _i  H  NO  4 

eg  eg  m  m  m  co  av  no  oo  on  n  oo  on 

•  •  . . 

oooo  ooo  ooo  ooo 


in  m  in 

m  oo  in 
n  <r  o 
in  non 


o 

in  o  oo 
m  m  o 
in  no  n 


N  fOOv 
no  eg  o 
vo  eg 
in  no  n 


00  O  O'  NO 
on  n  on  n 
in  in  on  no 
4  m  m  no 


ooo  ooo  oooo 


u  • 

pg  o 


3138  I 


O  -4  — s  -4  r--  Os 

2  — <  in  so  —i  in 

ooon^r^  co  n  n  cm  cn  cn 

cm  cm  cn  cn  so  —i  cm  as  — i  cm 

ooo . 


m  m  m 

— <  <r  so  o  — i  <r  so  H  4  so 


in 

as  — c 
cn  cm  — c 
C0  o  — c 


m  h  as 
masH 

co  o  N 

in  rv  rs 


in  in 

n  4  so  H  4  so 


OOO—i 


SO  CO  H 

so  so  cn 

os  4  CM 
so  oo  as 


in 

n  4  so 


oo  as  as 

O  4  CM 

cn  so  cn 
so  rv  oo 


m 

-i  4  so 
*  •  • 
O 


o  cm  cn  os 
in  oo  m  r-. 
4  m  — i  as 
in  so  n- 


in 

-i  4  X)  O 


OOO—l 


0.4931  0.4495  0.3389  1.289  0.8666  1.019 


-24- 


00 

00 

vO 

cn 

ON 

VO 

CN 

vO 

o 

ON 

in 

00 

CN 

00 

ON 

r- 

o 

o 

r- 

o 

o 

00 

in 

CN 

o 

r^. 

cn 

CN 

ON 

in 

cn 

o 

CM 

1  00 

o 

• 

On 

CN 

CM 

«— 4 

o 

o 

rH 

o 

o 

cn 

cn 

CN 

CM 

o 

i—4 

r-H 

i-H 

p—4 

CN 

rH 

CN 

*-H 

CN 

»-h 

CN 

00 

CM 

CN 

r- 

m 

cn 

nJ 

00 

o 

un 

cn 

cn 

CN 

in 

p—4 

o 

vO 

m 

ON 

3 

f-H 

o 

CN 

CN 

m 

NO 

cn 

r^. 

r- 

cn 

cn 

nO 

<r 

O 

O 

o 

ON 

NO 

o 

ON 

vO 

o 

H 

r- 

m 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

H 

o 

o 

o 

-o 

ON 

CM 

00 

-o 

•o 

in 

r*. 

in 

00 

CN 

in 

NO 

r— 4 

00 

NO 

ON 

CN 

o 

00 

in 

<?■ 

cn 

ON 

p—4 

00 

CN 

in 

cn 

o 

in 

r- 

r- 

On 

cn 

CN 

CN 

vO 

\o 

nO 

in 

cn 

cn 

CM 

CN 

-o 

cn 

r^. 

r^. 

NO 

r- * 

CN 

r-4 

H 

f— 4 

CN 

H 

p—4 

pM 

CN 

r-H 

r-H 

CN 

H 

CN 

00 

m 

00 

cn 

o 

<T 

00 

nO 

ON 

00 

00 

<t 

On 

r** 

ON 

00 

r-4 

p—4 

On 

ON 

00 

in 

cn 

On 

r^. 

MD 

00 

On 

»— 4 

CN 

f"- 

00 

ON 

r^. 

On 

o 

cn 

cn 

CM 

H 

o 

cn 

cn 

CM 

o 

cn 

CN 

rH 

o 

cn 

CN 

CN 

o 

cn 

CN 

o 

o 

o 

o 

o 

O  O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

NO 

m 

nO 

o 

CN 

r4 

iH 

ON 

cn 

cn 

H 

cn 

NO 

NO 

ON 

00 

cn 

cn 

o 

vO 

CM 

00 

r- 

cn 

CN 

nO 

cn 

NO 

r^. 

cn 

ON 

o 

00 

r^- 

cn 

00 

m 

00 

in 

ON 

m 

r- 

nO 

^4 

00 

NO 

m 

o 

VO 

CN 

ON 

cn 

ON 

in 

in 

vO 

cn 

<r 

m 

cn 

<r 

in 

in 

cn 

in 

m 

cn 

cn 

<r 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

cn 

vr 

r-4 

m 

ON 

m 

NO 

»—4 

ON 

CN 

o 

f—i 

cn 

NO 

m 

ON 

00 

r-- 

r*- 

vf 

CN 

m 

CN 

in 

CN 

-O 

in 

ON 

NO 

m 

un 

CN 

ON 

00 

nO 

o 

nO 

CN 

cn 

ON 

r^. 

CN 

ON 

p-H 

CN 

o 

o 

in 

vO 

r* 

cn 

<r 

in 

in 

<t 

m 

m 

NO 

cn 

<r 

un 

NO 

cn 

■o 

<r 

un 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

vO 

m 

m 

m 

in 

o 

<r 

NO 

o 

<r 

NO 

o 

H 

<r 

NO 

o 

H 

■o 

NO 

o 

o 

o 

o 

o 

o 

— 

o 

o 

o 

H 

o 

o 

o 

o 

o 

o 

f4 

-I  o 

O'  p—4  p—4 

00  r-)  l"- 

o  r-~  <r  o 


>J  ID  H  n 

ia  vj  n 

•  •  •  • 
H  H  H  CS| 


CM  m  CM 
00  On  00 


.1  ,1 


* 

U-I  o 


Table  2,  cont. 


I 


7230  0.7230  0.1965  1.220  0.4327  0.8654 

7979  0.7979  0.0  1.128  0.0  0.7979 


1 


Fig. 2  Variation  of  skin  friction  with  wall  temperature  for  large  0 
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Fig .3 — Variation  of  skin  friction  with  ts  for  large  >3 


i 


30 


Fig. 5— Variation  of  skin  friction  with  injection  for  large  0 
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